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I. INTRODUCTION 



The Standard Model (SM) gauge group is a priori consistent with a large flavor group 
that is not observed experimentally and Yukawa couplings and mixings are introduced to 
comply with such experimental evidence [I]. However, the pattern of fermion masses and 
mixing angles and the nature of the neutrino mass (whether Dirac or Majorana) is not 
constrained by the gauge symmetry and it is thus natural to go beyond the SM to seek for a 
theoretical explanation of the observed masses and angles and of their hierarchies. To date 
a promising direction towards the understanding of flavor physics in the SM - notably of 
fermion mass matrices and mixing - has been the study of particular textures (in the quark 
sector), such as the Nearest Neighbor Interactions (NNI) texture [2HU and the Fritzsch type 
textures [5] , and the introduction of flavor symmetry groups that constrain the structure of 
the fermion mass matrix to reproduce the desired texture. The mixing in the lepton sector 
has also given strong support to the idea of an underlying flavor symmetry and many flavor 
symmetry groups have been proposed so far [6 8]; of them discrete groups have proved to 
be quite successful and we keep exploring such direction in the present work. 

Although the aforementioned approach, involving (discrete) flavor groups, provides a pos- 
sible description of the pattern of fermion masses and mixing angles, it seems not to be able 
to explain the hierarchy (ies) among the different species of fermionic masses, i.e., the flavor 
symmetry explains the number and location of the "zeros" in the textures for the fermion 
mass matrices, but not the relative size among the non-zero entries. In other words it appears 
that, in order to obtain experimentally suitable values for such physical quantities, one needs 
to allow for parameters of quite different order of magnitude. Such hierarchy problems have 
proved to be quite hard, in the framework of four-dimensional non-supersymmetric theories. 
On the other hand in the context of extra-dimensions, several models - with flat or warped 
extra-dimensions - were proposed as solutions to hierarchy problems, mostly focusing on the 
gauge-gravity hierarchy. In particular the Randall-Sundrum model 1 (RSI) [9] has drawn 
a great deal of attention as it provides a possible description of the TeV-Planck hierarchy 
with a small (Planck size) extra-dimension and treats the TeV scale as a derived scale. The 
geometry of such space is a slice of five-dimensional Anti-de Sitter space, a warped product 
between four-dimensional space, and an extra dimension with the topology of an interval 
whose extrema are the fixed points of a Z 2 orbifold symmetry. The space thus "ends" on 
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two flat four-dimensional planes (branes) located at the extrema of the extra dimension. 
Extra-dimensional models, especially Randall-Sundrum models, have also been explored in 
the context of the SM physics as fermion masses might be explained in terms of the overlaps 
between the extra-dimensional profiles of fermions and Higgs(es) [TOHTl] . 

It is then interesting to consider the possibility of merging both scenarios - flavor sym- 
metries and hierarchy from extra dimensional settings - in a single setup and to determine 
if a viable model can be constructed. This work presents one such case. To do this we use 
a minimally-extended version of the SM that incorporates a discrete flavor group, and that 
reproduces fermion mass patterns as well as mixing angles in both the quark and lepton 
sectors. The model is minimal in the sense that it contains the SM matter fields only (in 
particular, it does not add right-handed neutrinos), an extended Higgs sector composed of 
two Higgs SU(2) doublets, a Lepton number violating scalar, and the smallest Abelian group 
that renders the NNI textures for the quark mass matrices, namely the cyclic group Z4. 

The model is detailed in the forthcoming section, however let us stress here that such 
approach was already taken by other groups. For example, in |15j . a realization of the 
so-called Lepton Minimal Flavor Violation scenario was obtained in the framework of RS 
with right-handed neutrinos. There it was found that the first Kaluza-Klein (KK) mass 
scale could be as low as 3 TeV. Along similar lines, but for the quark sector and providing 
a solution to the issue of flavor changing neutral currents in extra dimensional settings, the 
work in [16j extends the strong sector in the bulk with an additional SU(3), broken down to 
Quantum Chromodynamics (QCD) by boundary conditions that render the model safe, even 
for a KK scale of 2 TeV. Closer to our approach is the work by Kadosh and Pallante [T7] . 
where an flavor symmetry is introduced in a RS setup. The main difference with our 
model is that, unlike [T7j, we use a flavor group that reproduces the mixing matrices in both 
the quark and lepton sectors automatically, without the need of introducing additional flavon 
fields to generate non-zero entries and/or further suppressions. This makes our model more 
economical in terms of additional fields and more tractable, specially in the scalar sector. 
Another important difference is that in our model neutrino masses are generated radiatively 
and there is no need to introduce right-handed neutrinos, as mentioned above. 

The gauge and Yukawa sectors of our model are presented in section [XXJ where the ex- 
pressions for the mass matrices in both the lepton and quark sectors are described and 



explicitly shown. Section III shows the numerical analysis and results for both sectors and 
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some remarks about the scalar potential and flavor changing neutral currents. Finally we 
present the conclusions. An appendix has been included in order to maintain the discussion 
in the paper self contained. 

II. MODEL 

The five-dimensional field theory we consider is a higher-dimensional extension of the 
SM and lives in the RSI background. Such space is a slice of AdS 5 space of Planck-size 
curvature, where the fifth dimension y is taken to be an orbifolded circle Si/Z 2 of radius R, 
and fields are odd or even under the action of the orbifold, i.e. under reflection y — > —y. 
In other words the extra dimension has the topology of a line that stretches between y = 
and y = ttR and the latter are fixed points of the orbifold, where two flat 3-branes, one 
with positive and the other with negative tension, are accommodated. The brane located 
at y = is referred to as the "UV brane", whereas the brane located at y = ttR is referred 
to as the " IR brane" . The bulk line element is thus given by 

ds 2 = dy 2 + e-^ri^dx^dx" , -irR < y < ttR (1) 

with k ~ Mpi anck being (proportional to) the AdS curvature, and the TeV-Planck hierarchy 
is realized provided kR ~ 10. Such geometry is realized with a negative bulk cosmological 
constant and with the aforementioned brane tensions: all three are of Planck size. Moreover, 
all the fields we deal with have vacuum expectation values (vevs) of order at most TeV, and 
we may thus neglect their backreaction to the geometry Q. 

The five- dimensional model we consider has a gauge-symmetry group that includes the 
Standard Model custodial symmetry G C s = SU (3) c x SU (2) L x SU(2) R x U (1) b—l in order 
to suppress excessive contributions to the Peskin-Takeuchi T parameter |18. augmented 
with a discrete Abelian group Z 4 , hence G = Gcs x Z 4 . All the fields of our model are 
charged under Z4 and, with the exception of the Higgs fields, are all bulk fields, along the 
lines of what is done in [19]. This setup thus enjoys a very nice interpretation in terms 
of AdS/CFT correspondence as only the Higgs fields are seen as TeV-scale composites of 
the strongly-coupled sector of the four- dimensional theory and, since all other fields are 
bulk fields, all the phenomenology of the model is addressable from the weakly-coupled 
five-dimensional model. 



5 



The matter content of the model involves (three families of) left-handed quark doublets 
Qu and lepton doublets Lu along with two copies of right-handed quark doublets, Um = 



' URi 


and D Ri = 


(u, \ 


\di ) 




{ d Ri J 



and one copy of right-handed lepton doublets, E Ri 



\lmj 



the 



zero-modes of "tilded" fields will be projected out from the IR brane by Scherk-Schwarz 
gauge-symmetry breaking [19j and only their KK modes are eventually non-zero on the 
brane. We then use two Higgs fields $1 and $2 that are bi-doublets of SU(2)l x SU(2)r 
and are confined to the IR brane. Lastly, we introduce a bulk singlet scalar field h charged 
under hypercharge and lepton number that participates in neutrino mass generation. The 
Z4 charge assignment for all such fields is similar to the one given in [20J, namely 



(gi,52,g 3 ) = (2,0,3) 

(ui,u 2 , u 3 ) = (di, d 2 , d 3 ) = (3, 1, 2) 

for the quark sector, with = Q(Qu), Ui = Q(Um) and di = Q(Dm), 

(q;i, a 2 , a 3 ) = (2,0,3) 
(ei,e 2 ,e 3 ) = (3,1,2) 

for the lepton sector, with = Q{Lu) and = Q(Em), and 

(x) = i 



(2) 
(3) 



(4) 
(5) 



>1, <V2) 



1,2 , 



(6) 



for the scalars with fa = Q($j), and x = Q(h) = 1- The main difference with [20] is 
that the left-handed up-type and down-type quarks have the same Z4 charges; the reason is 
dictated by the fact that we are using Higgs bi-doublets confined on the IR brane. These 
fields transform as $' = e tXaLTaL $ e~ lXaRTaR under SU(2)l x SU(2)r gauge transformations 
(jaR. = T aL = \ a a-, where a a are the Pauli matrices) and, as explained later, are responsible 
for the spontaneous breaking of SU(2)l x SU(2)r — > SU(2)e> on the IR brane. In other 
words the vevs of are SU{2)r> singlets 



v 

T2 




1 



(7) 



so that = with $ = (ia 2 )$*(—ia 2 ), and Ur and D R thus enter on equal footing 

in the Yukawa lagrangian and must then have the same Z 4 charges. 
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It is easy to check that Q is invariant under the diagonal part of SU(2) L x SU(2) R , 
namely 

<5<$>= ? A Q (r i <$)-<$)r ii )=0. (8) 

and thus preserves SU{2) D . 

A. Gauge-symmetry breaking 



The bulk (custodial) gauge symmetry must be broken down to the SM gauge group on the 
four-dimensional effective action at TeV scale and usual electroweak spontaneous symmetry 
breaking must also be accounted for. We realize this in the same way as in [T7J [19] by 
orbifold Scherk-Schwarz projection, i.e. by assigning independent orbifold charges to a field 
at the two different ends of the interval, and with canonical spontaneous breaking. 

We can impose on the bulk fields a Z 2 x Z 2 charge assignment. For scalar fields this 
corresponds to the field transformations 



(f>(x, ~y) = Z(p(x,y) 
(x, tcR — y) = Z'(j)(x, 7tR + y) 



(9) 
(10) 



with Z,Z' = ± and y = 0,ttR being the fixed points of Z 2 x Z' 2 . For 5D spinors, single 
valuedness of the lagrangian upon the action of Z 2 x Z' 2 requires 



We can fix 75 



1 
-1 



ip{x,y) = Z"/ 5 ij;(x,-y) 

ip(x, nR + y) = Z'j B ip(x, irR-y) . 

and decompose the 5D 4-spinor as ip 



(ii) 



with £(rj) being 



left-handed (right-handed) Weyl spinors, in the 4D sense. At the fixed points, i.e. setting 
y — in (11), we thus have that the left-handed Weyl spinor £ has charges (Z, Z') and the 



corresponding right-handed Weyl spinor 77 has charges (-Z, —Z') [j] Below we mostly only 
care about the extra-dimensional zero-modes, whose profiles are summarized in Appendix [A| 



We thus have that the fermion left-handed zero mode (A12) only exists for charges (Z, Z') = 
(+, +) whereas the right-handed zero mode exists for (Z, Z') = (— , — ). In other words, other 

1 Henceforth, when referring to the fermionic Z2 x Z'^ charges, we will indicate the charges associated to 
the left-handed part of the doublet. 
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non-trivial orbifold projections lift the mass of the zero-mode. For vector fields, whose would- 
be zero-modes are constant, we again have that only (Z, Z') = (+, +) allows for massless zero 
modes and other charge assignments lift the mass and thus realize the low-energy symmetry 
breaking. 

In our model we break SU(2) R gauge symmetry down to U(l) R on the UV brane via 
Scherk-Schwarz mechanism, i.e. we assign (Z,Z') = (— ,+) to the electrically charged 
vector bosons of SU(2) R , and (Z, Z') = (+, +) to the neutral gauge boson. Since the former 
couple the upper and lower parts of the right-handed fermion doublets, single-valuedness 
of the bulk lagrangian requires that if the upper part is even, on the UV brane, the lower 
part must be odd and viceversa. Hence, in order to have zero-modes for both the upper 
and lower parts of the right-handed quarks, we need to double the number of right-handed 
fields [Tj5] in the quark sector, as already mentioned above. The charge assignment for such 
fields thus reads 



Qd+,+], u Rl = r i ' j , D Ri = r n J (12) 





for the quark sector, and 




L Li [+,+], E m = u ' (13) 



for the leptonic sector. A vev on the UV brane then provides the breaking U(1) R x 
U(1)b-l^U(1) y \M- 

On the IR brane the vevs of the Higgs bi-doublets provide the spontaneous breaking 
SU(2) L x SU(2) R x U(1) B -l -> SU(2) D x U(1) B -l- Hence finally the superposition of all 
such breakings only leaves U(l) em untouched. In fact 



(14) 



with r aD = r aL <S> 1r + 1l ® r aR being the generators of SU{2) D (above the tensor product 
is left implied.) 
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B. Yukawa Sector 



1. Quarks 



The mass terms for the quarks come from the following Yukawa interactions on the IR 
brane 

r=g S(y - irR) 



Yukawa 



/ 



dy 



A 



+ (llhQuix^y^x^Unj^y) + (^QUx^y^x^D^x^y) 



+ {i 2 d ) t3 Q Lt {x^y)^)D R] {x^y) 



+ /i.e., 



(15) 



where A = Mpi is the Planck-scale and (7„'^)ij are dimensionless parameters assumed of 
(9(1) that, together with the quark extradimensional profiles, generate the effective four- 
dimensional Yukawa couplings. The Z A charge assignments for the quark and scalar fields 
then induce the desired Yukawa NNI textures (writing separately each scalar contribution) 



^0*0^ 



* 
0* 



^0 ^ 



0* 
0*0 



(16) 



which after electroweak symmetry breaking lead to the following quark mass matrices 

/ n F« /ffi.\ n \ / n Td //f. \ n \ 



o ru®i) o 
Hfi<*i> o r« 3 ($ 2 ) 
^ o r« 2 ($ 2 ) ru®i) ) 



M r , 



o r d 12 {^) o 
r&^i) o r^ 3 ($ 2 
^ o r| 2 ($ 2 ) r-,($i> y 



(17) 



with denoting the vevs. The r"j d are the effective four-dimensional Yukawa couplings 
that depend on the fermion extradimensional profiles overlap with the Higgs bi-doublets at 
the IR boundary. Thus, under the zero mode approximation (ZMA) (see Appendix for the 
KK decomposition and the explicit fermion profiles), each of the Yukawa terms above looks 
like 



I-ttR 



dy e ^s(y - ^) 7tj ^(^(y) ^ (j) ^-(^-(y) c , fcfl 



A 



V2kR 



(1/2-4X1/2-4,) 



A^ [ e 2 ff fcfl(l/2-c«) _ 1 ][ e 2 W fcfl(l/2-c«y) _ X ] 



Q<3(x)*{x)Q<g>( x ), 
(18) 
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where q = U,D, $ = $1^, and i, j are family indices accordingly chosen (note also the 
inclusion of a canonically-normalizing factor for the Higgs bi-doublet). The effective Yukawa 
coupling, given by the product inside the curly brackets, depends on two quark localization 
c - parameters 



[lu,d) 



k 



(1/2 -4) (1/2 -4f) 



h3 A\ ^2^(1/2-4) _ ^2^(1/2-4'/) _ X ] 



(19) 



where c a Li = c d Li = since the left handed components of the u and d quarks form an 
SU(2)l doublet. In order to extract the c - parameters that lead to experimentally allowed 
observables, we follow the Harayama parametrization in reference [3J, which is a transfor- 
mation of the up - and down - type quark mass matrices to a basis such that these display 
the NNI form without modifying the mass eigenvalues nor the CKM matrix entries. Once 
parametrized in that form and following [20J, we make the assumption that the 1 — 2 and 
2 — 1 entries in M U)d are equal leading to 

^ q u ,d/y u ,d ^ 
M u4 = m ub q u4 /y u4 b u4 (y U)d ) 



\ 







d u ,d(yu,d) vl,d ) 



(20) 



where y Ujd are free parameters and 

m l,d + m l 

Pu,d = 



Qu,d 



m t,b 



t.b 



b Ul d(yu,d) 
d u ,d{y u ,d) 



I p u ,d + 1 - yt, d - Ru,d(yu,d) ( (h4 

^yu,d 



hu,d + 1 - ytd + Ru,d(yu,c 



Qu,d 

y u ,d 



(21) 
(22) 

(23) 

(24) 



with 



2 „.2 \l/2 



Ru,d(yu,d) = ((! + Pu,d - yt,d) 2 ~ A (Pu4 + qi, d ) + 8ql,dyl,d) 



M u ,d are real matrices arising from the phase factorization of M Uyd 



(25) 



M Uyd — P* d M u ^P Uyd , 



(26) 
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with P ud being diagonal phase matrices such that P = P u P d = diag(l, e l ^ ud , e taud ), (5 u d = 
j3 u — (3d and a u d = a u — ad- Therefore, four parameters y U)d , Pud and a u d have to be chosen 
to fit the CKM matrix, 

Vckm = OlPOd, (27) 

where Ut d diagonalize M u>d M^ d 

Ol d M u , d Ml d O u4 = diag(m* , m c 2 m 2 ). (28) 



We show the numerical results that reproduce the experimental values in section III 



2. Charged Leptons 



The charged lepton masses are similarly obtained from the Yukawa interactions on the 
IR brane 



Yukawa 



dy- 



-g 8(y- % R) 
A 



(^Luix^y^ix^Enjix^y) 



+ ( 7e %.Z Li (^, y)<S> 2 (x»)E Rj (x>*, y) 



(29) 



After EWSB we obtain 



( rfa^i) 

V ^($2) r| 3 ($i) ) 



(30) 



The analog of Eq. (19) for charged leptons is 

k 



(7.) 



1/2 -eg,) (1/2 -tf 



'■j 



(31) 



'A y [ e 2 7 rfci?(l/2-c£ i ) _ ^2^(1/2-^.) _ ^ 

We can parametrize the M e matrix ^ following again [3], and so M e has dependence on both 
dimensionless parameters y e and z e 



( 



q e z e /y e o 

q e / (y e z e ) y/B e (y e ,z e ) 



\ 



y y/D e (y e ,z e ) 



lit 



(32) 



2 We do not assume (M e ) 12 = (M e ) 2 i as we did in the quark sector. 
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,2^2 



where p e and q e are the analogs of p u ^ and q u ^ for charged leptons and 

B e (y e , z e ) = \ {p e + 1 - yj ± y/(l -p e + yty - %| - y*z*)(q* - y*/z*)} - 

D e (y e ,z e ) = \ \p e + 1 - yt T ~Pe + yi) 2 ~ We ~ VMM ~ vl/4)} - ~ -(33) 
2 1 J y l e 

Observe we can have two different solutions: the plus (minus) case is obtained taking the 
plus (minus) and minus (plus) signs in B e and D e respectively. Depending on which solution 
we choose we have two different regions for the mathematically allowed values of z e and y e , 
Region I (plus case) and II (minus case) respectively. Each region is constrained by the non 
- negative real values of B e an D e . 

The matrix M e allows us to obtain the lepton c - parameters by comparison to Eq. (30). 
Under the assumption that the phases in the charged lepton sector are zero, there are only 
two parameters left, y e and z e , whose actual values will be set by the neutrino sector results. 



3. Neutrino sector 



Since there are no right-handed neutrinos present in the model, neutrino masses are 
generated radiatively ^ as in [20]. In the present extradimensional setting this mechanism 
is going to set two of the lepton c^-parameters. As can be noted in the next section, their 
values turn out to be 0(1), as any localization c - parameter must be. 

The radiative mechanism introduces a cubic Lepton number violating scalar interaction 
among the two SU(2) l doublets and the charged scalar h, as well as the Zee operator coupling 
the left-handed lepton doublet to the singlet charged scalar [2HI22]. r ^^ ie R an dall-Sundrum 
geometry enhances this cubic scalar interaction as well as the Zee operator with factors 
proportional to the extradimensional profiles of the fermions and, if allowed to propagate 
through the bulk, those of the scalars. In the scenario where both scalar Higgs doublets are 
confined to the IR brane and the scalar singlet h is a bulk field, the cubic scalar operator is 
an interaction on the TeV brane which looks like 

- Cmh = J dyy/^g~5(y - Tri^A^e^^)^ (x)h(x, y) , (34) 

3 Using the well known dimension five Weinberg operator leads to unacceptable large values for the neutrino 
masses. This could be remedied by allowing unnatural small values of the dimensionless 7 parameters, 
in contradiction to the philosophy of the general scenario. We therefore stick to the 4D renormalizable 
argument in |20) and consider the radiative mechanism. 
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where i, j are SU(2) indices and \ a p is antisymmetric and with mass dimension +1/2. Under 
the ZMA approach this term acquires the form 

- = r R dy e-^5(y - 7riZ)A a ^€ -^(g)^(g) fe( ° )( ^ ^(e' fcfl ) 2 , (35) 

with the dimensionless antisymmetric coupling Aq,^ of 0(1) (note also the inclusion of a 
canonically-normalizing factor for each one of the Higgs doublets). Thus, for an UV-peaked 
bulk h field (whose approximated profile and normalization factor are shown in appendix [A| 
the effective cubic operator acquires the form 

- C** h = X a ^t(x)^(x) ■ r]f V) , (36) 

where 



{uv) - *fky/k\b. - l\e^ b - 2 ^ kR (37) 



has the correct mass dimension: [t^ 7 ^] = +1 and 6_ < 1. An UV-peaked bulk h (again in 
the ZMA) also modifies the Zee operator 

C LLh = J dy^k ab e lJ (Ll i nx,y)L b L3 (x,y)h*(x,y) , (38) 
(where [k ab ] = — 1/2 and antisymmetric by the Pauli principle) and leaves it as 

C LLh = I dy e-^lvl « eg(Lff> )c {x)L m { x )^> {x ) l^M UjW h^M , (39) 



tR 



VX V " Li ' K ' Lj " J " J V^R V^R V^R 



where K ab is now dimensionless (and antisymmetric). The effective Zee operator then looks 
like 



- C LLh = ^e l3 {LT)^)L b ^{x)h^\x) ■ ^ , (40) 
where r]^^ is dimensionless and given by 



AUV) = k ^lJ (5 C La)(\ C Lb)\ h - l \ f +7TR , (6--cf; o -c£ 6 )A[«| (m 



J W [ e 2rffi(l/2- C iJ _ 1 ][ e 2 7 rfeii(l/2- C ^) _ !] 



Following the expressions in [20] for the 4D Majorana neutrino mass matrix and replacing 
the effective Yukawa and scalar couplings, we find that the Majorana neutrino mass matrix 
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entries are 

m VeVe = 2a e K 31 V ^m T ^\ 12 r ] £ rv \$ 2 )F(ml,m 2 h ), 
= 0, 



m 



(42) 
(43) 
(44) 

m v ^ = m^ Ue = 2b e K^ v) m TT [\ 21 r)V v \<b l )F(ml 1 m\), (45) 
m VeVr = m VTVe = 2^a e K 13 V ^m ei ,\i2V { x V) ^2) + b' e K 31 rj K m T ^X 21 rj^ V) <<£>!> (46) 

+ c e K^ v) m TT \i2vf V \^F{ml,ml), (47) 
0, (48) 



to 



TO,. 



where F(to|,to^) (accounting for the scalar loop factor) is a function depending on the 
scalar masses of the charged Higgs, to$, and of the singlet scalar m h through 



F(ml,m 2 h ) = 



1 



log 



tos 



16n 2 to| — to| ml 



(49) 



The parameters a e , b e , b' e , and c e belong to the Yukawa matrix for charged leptons 



/ 



\ 



V\ 2 

1 21 u 1 23 
u 1 32 1 33 



\ 



a e 
al b P 



\ 



&' e C e 



\ 



and the entries to^, i, j — e, /i,t are those of the matrix M e . 

Now, since f% and Ay are antisymmetric we can rewrite the neutrino mass matrix entries 

as 



m UeUe = - (m efl m TII tan (3)C, 
m v u v u = 0, 



tan (3 



tan (3 



m UeUT = m UTVe = ^TO^tan/3 + j^-^ - m^tan/3 ) C- 
m ^u T = = 0, 

where C is a common factor with dimension of inverse of mass 

A3^UV\ JJV 



C = 2K l X V *i2V U x V F(rnl,rnl), 



(50) 
(51) 
(52) 

(53) 

(54) 
(55) 

(56) 
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and 



tan (5 



(gg) 



(57) 



III. NUMERICAL ANALYSIS 



A. Lepton sector 



To perform the numerical analysis in the lepton sector we used the experimental data at 
3<7 from the global neutrino data analysis in 





Best Fit Value 


3(T range 


sin 2 9i2 


0.312 


0.27-0.36 


sin 2 6*23 


0.52 


0.39 - 0.64 


AmlJlO^eV 2 ] 


7.59 


7.09-8.19 


Am 2 2 [10~ 3 eV 2 ] 


2.50 


2.14-2.76 




-(2.40) 


-(2.13-2.67) 



(58) 



with 5c p = and normal (inverted) hierarchy and the recently Daya Bay results (confirmed 
at 5a) [21] 

sin 2 29 13 = 0.092 ±0.017, (59) 

which can be rewritten as 

sin 2 #i3 = 0.0235 ± 0.0045. (60) 

By convenience in the analysis we also define the following range for the neutrino squared 
mass differences ratio (at 3cr) 



NH (IH): 0.027(0.028) < 



Am 2 ! 



Am 2 2 



< 0.035(0.036), 



(61) 



obtained by summing in quadrature the relative errors of Am^i and |Am 2 2 |. The con- 
straint over the sum of neutrino masses, ^2m Vi < 0.29 eV presented recently in [25] is also 
considered. 

The charged leptons masses used in the analysis are those given by the central values 
in [26] 



m P = 0.511 x 10" 3 GeV, 



m M = 0.1056 GeV, 



m T = 1.776 GeV. 



(62) 
(63) 
(64) 
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As in reference [T7] we take the 5D scale to be A = k = M P i anck where M P i anck = 
2.44 x 10 15 TeV is the reduced Planck mass, and the effective scale ke~ nhR ~ 6.89 TeV 
generated by irkR = 33.5. For the other parameters involved in the neutrino mass sector we 
take m$ = 0.5 TeV (recall this is the mass of the charged Higgses in the loop), rrih = Mpi anc k, 



K 



1 ' — 1 = — k 31 , and A 12 = 1 = —A 21 . Note that since h is UV-peaked its mass is taken to 



its natural value of Mpi anc k. We note however that the model can reproduce the neutrino 
sector even in the case of a very light /i-field (or even an IR-peaked h field) due to the fact 



that to/i only enters through the common factor C in Eq. (56), and thus does not affect the 
diagonalization. Its only possible effect is in the absolute size of the neutrino mass matrix 
entries and it corresponds to acceptable changes of 0(1). 

From the expressions for the neutrino mass matrix entries we observe that the only 
coefficients involved are and c£ 3 , and are taken to be c\ x = 0.8 and c£ 3 = 0.55 (this is 
only a choice and corresponds to similar values used in [17]). We also observe that all these 
parameters are contained in the C factor, together with the parameter b = 2—^4 + a_, and 
in consequence they do not affect the diagonalization matrix. Thus, the lepton mixing matrix 
Upmns defined by Upmns — U\U V , where Up and U v are the diagonalization matrices in 
the charged lepton and neutrino sectors respectively, only depends on tan (3 and the charged 
leptons mass matrix entries. 

We use the standard parametrization of the Upmns given in [26] 

V = U PMNS P, (65) 

where 



UpMNS 



( C13C12 C13S12 s 13 e l5cp ^ 

-C23S12 - S 2 3S 13 C 12 e lScp C23C12 - S 2 3S 13 S 12 e lScp S23C13 
y S23S12 — c 23 s 13 c 12e* <5c ' P — S23C12 — C 2 3S 13 Si2e*' 5cp C23C13 J 



(66) 



with Cij = cos 8ij, Sij = sin 6^ and where Sep is the Dirac-CP violating phase in the neutrino 
sector which we assume is zero. P = diag(l, e ta21 , e*" 31 ) is a diagonal phase matrix containing 
the two Majorana CP violating phases. 

The first step in our analysis consists of a scan over the mathematically allowed regions 
for the z e and y e parameters of the charged lepton mass matrix M e that lead to positive 



values for B e and D e in Eq. (33). The next step is to scan over values for tan/3 values that 
lead to the allowed neutrino mass differences ratio and mixing angles. Finally, once the 
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ranges for z e , y e and tan f3 values that give the right angles and ratios are determined, we fit 
the a_ values (or b_) required to obtain the neutrino masses satisfying the constraint over 
the sum of neutrino masses [25] . 

Performing the scan over tan/? values from 0.1 to 1.0 (in steps of 0.1) and from 1 to 20 
(in steps of 1) we find that in Region I (plus case) the only value that works is tan = 2.0 
while in Region II (minus case) tan/3 = 0.3 and tan/3 = 0.6 give acceptable values. 

Figure [T] shows the parameter space obtained in this process. Note that in all cases there 
is a large region in the z e — y e plane consistent with positive values for B e and D e . However 
the allowed region, consistent with all experimental information, is considerably reduced to 
small ranges in this parameter space. In the plots we label by allowed region the region of 
parameter space consistent with all experimental data, including the constraint on the sum 
of neutrino masses [25J. 

We note that, as found in [20], only inverted hierarchy for the neutrino masses is obtained 
in this model and that the resulting Majorana phases are found to be 0121 = «3i = tt '. 

It is important to note that throughout the analysis all the 7?- have been taken to be of 
(9(1). In fact, except for jf 2 and 7f 2 which are used in the fit and take values of 0(1) in all 
cases, all other are set equal to 1. The same is true for all the profile c - coefficients: all are 
of 0(1). 



B. Quark sector 



In [27] the values for the parameters y u ^, (3 u d and a u d leading to a correct CKM matrix 
\Vckm\ were found to be y u = 0.9964, yd = 0.9623, (3 u d = 1.4675 and a U d = 1.9560. The 
corresponding mass matrices (in TeV) read 

/ n n nnnn^^ns n \ 



M„ 



0.0000552118 

0.0000552118 0.0118527 

0.0183451 0.17152 



( 



0.0000235694 

0.0000235694 0.000268666 

0.00158321 0.0039362 



(67) 



(68) 



Comparing these matrices with Eq. (17) and using the values for tan (3 found to work in the 



neutrino sector, we performed a fit to the c - profile parameters and 7's. We find that it 
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Figure 1. Consistent regions of parameter space in the z e — y e plane for Case I and tan/3 = 2 
(upper plot), Case II and tan/3 = 0.3 (middle), and Case II and tan/3 = 0.6 (bottom). We show 
the complete region consistent with positive values for D e and B e , the region consistent with the 
neutrino mass squared differences ratio, the regions consistent with neutrino mixing angles, and 
the allowed region (see text). Note that for Case II and tan/3 = 0.3 the region consistent with 
mixing angles is basically the same as the allowed region. 
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is possible to obtain solutions with all parameters of 0(1). In particular, for tan/3 = 0.3, 
experimental agreement is found for at least the following set of values (c^, C£ 2 , c^) = 
(0.7,0.491253,0.128276), {c u R11 c u R2 , c u m ) = (0.666249,0.342679,0.737729), (c d R1 , c d m , c d m ) = 
(0.690378,0.468589,0.591601), (7H) = l 3x3 , 7^ = 0.101081, 7^3 = 2.02488, and all other 7's 
equal to 1. 

Thus, using dimensionless parameters of order one and dimensionful parameters associ- 
ated to the scale of the setup, we find that it is possible to reproduce all observed masses and 
mixing angles in both the quark and lepton sector. Furthermore we accomplish this with a 
minimal set of additions and only left-handed neutrinos. It is important to recall that all re- 
sults presented in this work have been obtained under the ZMA scheme and thus neglecting 
(although safely [TTJ ) possible contributions from KK modes. Yet another possible source 
of contributions that must be investigated corresponds to contributions from higher order 
operators consistent with the gauge and flavor symmetries. In our case, the smallest higher 
order operators that could contribute are of the general form $$$QlQ,r (and similarly for 
leptons). We find that the gauge and flavor invariant contributions, denoted by Sm^, to the 
zero entries in the mass matrices satisfy Smij/mik ~ 10~ 7 — 10 -11 , where the denote the 
values of the non-zero entries. For the case ij = Ik, corresponding to contributions to the 
non-zero entries, we find 5rriij/mij ~ 10 -4 . Thus, these contributions can be safely ignored. 

Before concluding we make some remarks to be taken in full consideration in a future 
work. Our model has a potentially interesting scalar phenomenology that requires a complete 
study, including its vacuum stability and possible collider signals. Another important and 
related phenomenological issue is the presence of flavor changing neutral current (FCNC) 
effects. In our model there are two different possible contributions to FCNC that must 
be taken into consideration. The first one comes from the fact that we are dealing with a 
two (flavored) Higgs doublet model in which each Higgs couples to both the up and down 
quark sector. This can result in tree-level FCNC [2HI EH] that need to be analyzed and will 
induce constraints on the parameter space in the scalar potential. The second possibility 
is associated to the presence of mixing between the degrees of freedom of the effective 4D 
theory and their KK excitations. This is a problem that all extra dimensions models have 
to face. One has to guarantee that the corrections induced by this mixing will keep the 
agreement with observations forcing a lower bound, of order a few TeV, on the KK scale. 
In our model, the use of a custodial symmetry allows to consider a lower bound on the first 
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KK mass to that imposed by electroweak precision measurements [TS] (see also [T51 EES1 EE] ) ■ 
Furthermore if the flavor pattern of the Yukawa couplings holds in the 5D theory due to 
a flavor symmetry, an increased alignment between the 4D fermion mass matrix and the 
Yukawa and gauge couplings is obtained suppressing the amount of flavor violation induced 
by the interactions with KK states. Moreover, as pointed out in [17], a bulk flavor symmetry 
might also induce a cancellation of observable phases, therefore the dominant new physics 
contributions to the neutron and electron dipole moments and/or to e'/e# will vanish. 



IV. CONCLUSIONS 



A five- dimensional warped extension of the Standard Model is considered, where the 
space-time background is taken to be the Randall- Sundrum model 1, i.e. a warped extra- 
dimension of Planck size in an Anti-de Sitter bulk. We consider a Zero-Mode Approximation 
for the four-dimensional effective action one obtains as a dimensional reduction of the higher- 
dimensional theory and that, in the present case, constitutes an extension of the SM. In 
the model we study above almost all the fields in the 4D effective action are zero modes 
of bulk fields, in particular each left- respectively right-handed fermion mode has its five- 
dimensional counterpart; only the Higgs fields are purely four-dimensional and live on the IR 
brane (i.e. on the 4D plane located at the orbifold fixed point y = ttR) and, unlike the bulk 
fields, have no KK excitations. Bulk fermion fields come about with a (Planck size) mass 
term parametrized by a real number c which determines the shape of the extra-dimensional 
profile of the field. The Yukawa interaction between Higgses and fermions is given in terms 
of overlapping integrals of zero-mode profiles for the fields involved and thus depends upon 
the aforementioned c real parameters along with numerical coupling constant (hermitian) 
matrices 7's. In order to obtain a NNI-type quark mass matrix, we advocate the presence 
of a discrete symmetry upon which all fields are charged: the cyclic group Z4 is the smallest 
Abelian group consistent with the aforementioned texture. Also for the charged lepton sector 
the discrete symmetry fixes the form of mass matrix in a similar way as for the quark sector. 
For the quark sector (and similarly for the charged lepton sector) the numerical analysis goes 
as follows: first we consider the mass matrix parametrized d la Harayama in order for it to 
display NNI form. The diagonalization of the real squared counterpart of the previous mass 
matrix leaves four parameters that can be fixed using the experimental values for the CKM. 
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In turn, this allows to fix the entries of the mass matrix itself and, by comparison of such 
entries with the ones foreseen by the overlaps of extra-dimensional profiles, one can obtain 
the compactification parameters c's and the Yukawa dimensionless entries 7's. All these 
parameters turn out to be of 0(1). For the charged leptons the parametrization goes slightly 
different and the (real) mass matrix depends upon two free parameters only constrained by 
the reality of the mass matrix: allowed parameter regions satisfying such constraints, as well 
as all experimental data, are graphically displayed. In the neutrino sector we choose not 
to introduce right-handed neutrinos in the model and generate neutrino masses radiatively. 
This is accomplished with the help of a bulk charged scalar field whose zero-mode is peaked 
on the UV brane (located at the fixed point y — 0): such field also mediates the lepton 
number violation. With scalar field masses chosen to be m$ = 0.5 TeV for the charged 
Higgses in the loops, and rrih = M P i anck TeV, experimental data taken from 3a global 
neutrino data analysis and Daya Bay results are matched by setting all c parameters and 7 
parameters to values of 0(1). Although a complete phenomenological study of the model 
is under preparation, some comments regarding the scalar phenomenology of the model as 
well as its possible contributions to FCNC were also briefly discussed. 
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Appendix A: On the zero mode profiles 

We briefly review the computation of zero- mode profiles for some bulk fields, on the RSI 
background considered above; we basically follow what is done in [13], [28]. A bulk field 
satisfies a second order differential equation of the form 



e^rTd^du + e sa d y (e~ sa d y ) - M, 



$(x,y) = , 



(Al) 



where $ = 0, e- 2a ^ LtR }, s = {2,4,1} and M| = {0,ak 2 + ba'\c(c+ l)k 2 - ca"}; 
note that compared to p31 [28] 's we redefined c — > — c for the right-handed mode. Above 
c(\y\) = k\y\ (modulo 2tiR) and thus o' = ke(y) and a" = 2k[5(y) — S(y — irR)]. Upon 
Kaluza-Klein decomposition one gets 



1 oo 



(A2) 



n=0 



with rj^d^d^^ix) = m^^(x), and orthonormality conditions 

e 2lsr7 f n (y)fn>(y) = 5 ny , Z s = (0, -1, -3/2) 



27ri? 



so that 



fn(y) = m 2 J n (y) . 



(A3) 



(A4) 



Here we only concentrate on the zero modes for which mo = 0. 

Vector field 

For the vector field it is immediate to realize that the only possible zero-mode profile is 
constant and since the measure is also trivial in this case, we simply have fo(y) = 1. 

Scalar field 



For a bulk scalar field we have 



-e Aa d y {e- ia d y ) + ak 2 + ha" 



fo(y) = o 



that admits solutions 



1 



fo,±(v) = ^ b±a > K = 2 ± , 

iV 



(A5) 



(A6) 



22 



and the normalization factor is given by 



N , 



e 2(b ± -l)-wkR 
27TkR(b ± - 



1 1/2 



1 



3 (6-l)7TfcR 



b = b + 
b = b_ 



(A7) 



'2irkR\b-l\ i 

where in the last expression we took into account that krcR » 1 and b + > 1, b_ < 1. 
Hence, including the measure factor the zero-mode profiles read 

,(6 + -i)fc(|v|-7rJi) b = b. 



1 



- a fo(y) « v4&^ 



,(6_-l)fc|»l 



(A8) 



so that for b = b + the profile is peaked about the IR brane (y = ttR) and for b = b_ the 
profile is peaked about the UV brane (y — 0). 



Fermionic field 



For a bulk fermionic field we follow [13] with the aforementioned renaming of parameter 
c for the right-handed mode. Hence 



-e a d y (e~ a d y ) +c(c+l)k 2 - ca" 



e~ la k{y) = 



and 



with 



foL,n(y) = TT e 

iVn 



(2-c)a 



Nn 



e 2nkR(l/2-c) _ l^ 1 / 2 

2nkR(l/2 - c) 



1 



e (l/2-c)irkR ; c< i/2 

c > 1/2 



'2<KkR\l/2- c\ i j 

Hence, including the measure factor the zero-mode profiles read 

,(l/2-c)k(\y\-7rR) ^ c < 1/2 

c > 1/2 



1 



,-3o 



. ^f (y) » >/|l/2-c|fc 

V27T.R | g(l/2- C )%| 



(A9) 



(A10) 



(All) 



(A12) 



so that for c < 1/2 the profile is peaked about the IR brane (y = tcR) whereas for c > 1/2 
the profile is peaked about the UV brane (y — 0) . 
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